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1.0  INTRODUCTION 

The  use  of  laser  scattering  techniques  such  as  spontaneous  and  nonlinear  Raman 
scattering  for  diagnostics  of  gas  flow-field  and  combustion  environments  has  been 
successfully  demonstrated  previously  by  various  groups.  This  particular  panoply  of 
techniques  provides  the  capability  to  perform  local,  time-resolved  measurements  of  the 
number  density  of  individual  species  as  well  as  the  energy  distribution  or  enthalpy  of  the 
rotational  and  the  individual  vibrational  modes  of  molecular  species.  The  previous 
applications  of  these  measurement  methods  have  included  measurements  of  the  free-stream 
flow  properties  of  nozzles  (Refs.  1  and  2)  and  orifice  sources  (Ref.  3)  as  well  as  studies  of  the 
boundary  and  shock  layers  of  an  aerodynamic  model  (Refs.  4  and  5),  measurements  of  the 
exit  plane  properties  of  monopropellant  liquid  rocket  exhaust  plumes  (Ref.  6)  and  a 
turbulent,  two-stream  H2-air  combustion  flow  field  (Ref.  7),  and,  more  recently, 
spontaneous  and  nonlinear  Raman  scattering  studies  of  both  spark  ignition  engines  (Refs.  8 
and  9)  and  compression  ignition  engines  (Ref.  10).  The  success  of  these  diverse  applications 
has  been  such  that  the  focus  of  research  and  development  in  several  areas  has  shifted  from 
the  question  of  feasibility  of  the  techniques  to  that  of  the  ultimate,  but  practical,  uncertainty 
of  the  methods  of  measurement.  To  appreciate  the  question  of  ultimate  uncertainty  of  the 
measurement,  it  is  helpful  to  recall  the  basic  equation  which  relates  the  scattered  signal  S  to 
the  power  P0  of  the  laser  of  frequency  t>0,  the  subtended  collection  optics  solid  angle  (AQ), 
the  total  number  density  (n-r)  of  the  scattering  species  and  the  thermally-averaged 
differential  cross  section  (do/dli)  for  the  transition  of  interest.  Although  the  signal-to-noise 
ratio  is  the  pertinent  parameter,  for  the  purpose  of  this  work  it  is  sufficient  to  note  that 

S  =  nT  )  x  x  Afl  x  t 

'  dO  *  v  h*0  / 

from  which  it  is  seen  that  S  varies  directly  with  the  differential  cross  section  do/dfi  . 
Therefore,  knowledge  of  dtx/dO  for  the  species  is  essential  for  quantitative  determinations 
of  density  n  and  rotational  and  vibrational  temperatures  TR  and  Tv,  respectively. 

If  the  gas  sample  of  interest  can  be  provided  for  a  laboratory  calibration  for  which  the 
conditions  can  be  varied  to  span  the  range  of  parameters  of  the  unknown  gas  sample,  it  is  a 
simple  matter  to  perform  an  in  situ  calibration  and  to  use  theoretical  predictions,  which  have 
been  verified  by  the  laboratory  calibrations,  to  determine  the  desired  properties  of  the  gas. 
However,  it  is  not  uncommon  that  the  properties  of  the  unknown  gas  sample  cannot  be 
spanned  by  conventional  laboratory  studies,  and,  consequently,  the  theoretical  predictions 
are  unverified.  For  simple  and  well-understood  species  such  as  N2  and  H2  such  an 
extrapolation  is  normally  justifiable.  However,  with  increased  complexity  of  the  molecular 
species  arise  difficulties  in  the  interpretation  of  spectral  intensities.  These  difficulties  are 
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typified  by  complex  overtone  spectra  and  intramolecular  perturbations  such  as  the  Fermi 
resonance  phenomenon.  For  such  molecular  species,  at  the  high  vibrational  temperatures 
encountered  in  various  applications  of  spontaneous  and  coherent  Raman  scattering,  these 
uncertainties  in  knowledge  of  the  spectral  intensity  distributions  are  reflected  as 
uncertainties  in  the  Raman  cross  sections  and  result  in  uncertainties  in  the  measured  values 
of  species  density  and  temperatures.  Therefore,  calibrations  or  cross  section  measurements 
are  required  for  these  species  over  the  temperature  range  of  interest  for  the  required 
applications. 

The  temperature  range  of  interest  for  the  various  applications  of  Raman  scattering  is 
quite  broad  and  spans  the  nominal  range  of  20-3000  K.  Calibration  and  cross  section 
measurements  for  the  20  to  =  1200  K  range  can  be  performed  using  a  combination  of 
experimental  methods  including  expansion  flow  fields,  cryogenically  cooled  gases, 
resistance-heated  furnaces,  and  compression  ignition  engines  which  are  operated  in  a  non¬ 
combusting  mode  (Ref.  10).  However,  for  the  higher  temperature  range  (1200-3000  K)  the 
choices  are  less  abundant.  Even  though  a  variety  of  techniques  exists  to  heat  gases,  for  cross 
section  measurements  it  is  essential  that  the  gas  be  in  a  state  of  local  thermal  equilibrium 
(LTE),  at  least  with  respect  to  the  modes  of  molecular  motion  involved  in  the  Raman 
process.  Because  of  this  requirement,  the  available  methods  for  providing  the  LTE  gas 
sample  are  reduced  to  essentially  two:  the  flat-flame  burner  (Ref.  1 1)  and  the  shock  tube 
(Ref.  12).  The  burner  method  has  been  used  successfully  by  Lapp  and  co-workers  (Ref.  13) 
for  a  variety  of  cross  section  measurements  of  gases  which  occur  normally  as  a  result  of  the 
combustion  process.  However,  there  are  features  of  this  approach  which  contraindicate  its 
use  for  general  cross  section  measurements.  First,  the  temperature  range  of  the  flat-flame 
burner  is  inadequate  because  3000  K  cannot  be  obtained  routinely.  Second,  unless  the 
species  of  interest  is  a  normal  combustion  product  of  the  flame  reactants,  the  selection  of 
which  determines  the  flame  temperature,  the  species  must  be  injected  into  the  flame  with  a 
concentration  sufficient  for  the  cross  section  measurement.  However,  this  very  injection  can 
alter  significantly  the  properties  of  the  flame,  thereby  invalidating  its  status  as  a  standard. 
Additional  flame  calibrations  will  then  be  required  for  each  concentration  level  of  the 
injected  species,  and  the  precision  of  the  cross  section  measurement  is  degraded.  Moreover, 
it  is  possible  that  the  injected  species  will  experience  chemical  reactions  with  the  flame 
species,  further  complicating  the  task.  Finally,  for  a  burner  of  reasonable  size,  the 
magnitude  of  the  spatial  gradients  of  the  flame  properties  is  a  serious  consideration,  and, 
quite  likely,  these  gradients  result  in  severe  restrictions  in  the  experimental  design  of  the 
cross  section  measurements. 

The  shock  tube  alternative  to  the  use  of  combustion  flame  gas  sources  is  a  relatively 
simple  laboratory  tool  for  the  production  of  high-temperature,  LTE  gas  samples;  a  sketch  of 
a  typical  shock  tube  is  shown  in  Fig.  1.  For  ideal  operation,  the  desired  high-temperature  gas 
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sample  can  be  provided  over  a  range  of  temperatures  by  simply  varying  the  ratio  of  the 
initial  driver-driven  gas  pressures  (Ref.  12).  Also,  the  shocked  gas  density  can  be  varied 
while  the  shocked  gas  temperature  is  held  constant  by  simultaneously  varying  the  driven  gas 
pressure  and  the  ratio  of  the  driver-driven  gas  pressures.  Flexibility  in  the  experimental 
design  is  achieved  by  observation  of  the  even  greater  temperature  and  gas  density  values  of 
the  reflected  shock  region.  Additional  advantages  of  the  technique  include  the  capability  of 
operation  of  the  shock  tube  with  the  species  of  interest  as  a  pure  gas  sample,  unlike  the 
combustion  source,  and,  using  a  reasonable  diameter  shock  tube,  measurements  can  be 
performed,  in  principle,  with  an  LTE  gas  sample  with  small  spatial  gradients.  Therefore,  as 
will  be  shown  in  the  following  section,  the  shock  tube  provides  a  method  of  spanning  both 
the  temperature  and  number  density  ranges  required  for  the  various  cross  section  studies. 

The  major  disadvantage  of  the  use  of  shock  tubes  for  cross  section  measurements  has 
been  the  comparatively  short  run  times  and  data  acquisition  periods  available  for  the 
measurement.  For  example,  10-50  ftsec  represents  a  typical  interval  of  time  available  for- 
measurements  in  the  incident  shock  region;  operation  of  the  shock  tube  in  the  reflected 
mode  can  increase  this  interval  by  over  an  order  of  magnitude  (Ref.  12).  Regardless  of  the 
mode  of  operation  of  the  shock  tube,  data  acquisition  times  of  50  jisec  -  1  msec  are  quite 
short  in  comparison  to  the  time  available  for  continuous  combustion  sources.  Clearly,  the 
availability  of  pulsed  laser  sources  lessens  the  severity  of  the  short  run-time  disadvantage, 
and  the  use  of  a  Q-switched  laser  source  with  a  nominal  10-nsec  pulse  width  easily  achieves 
the  required  temporal  resolution.  However,  the  flow  duration  disadvantage  is  not  eliminated 
because  of  the  methods  normally  employed  for  spectral  dispersion  and  subsequent 
detection.  Specifically,  detection  and  registration  are  desired  for  the  spectrally  resolved 
scattered  radiation  for  broad  wavelength  regions,  but  practical  requirements  for  the 
minimum  acceptable  signal-to-noise  ratio  (S/N)  of  Raman  scattering  and  the  typical  pulse 
duration  and  repetition  rate  of  the  laser  source  prevent  the  use  of  rapid-scan  spectrometers. 
The  utilization  of  single  or  multiple  photomultiplier  tubes  (PMT)  with  a  fixed  wavelength 
setting  of  the  spectrometer  is,  of  course,  possible  (Ref.  10).  However,  using  this  approach  to 
cover  the  desired  broad  spectral  region  and  to  simultaneously  provide  the  spectral  resolution 
necessary  for  the  cross  section  measurements  requires  the  use  of  multiple  shock  tube-laser 
source  firings  in  conjunction  with  variation  in  the  central  wavelength  of  the 
spectrometer.  Consequently,  because  of  the  shot-to-shot  variations  in  both  laser  and  shock 
tube  parameters,  the  resulting  cross  section  data  can  be  expected  to  exhibit  larger  than 
necessary  experimental  uncertainty  values.  As  a  result  of  these  anticipated  difficulties,  the 
use  of  shock  tube  flows  for  Raman  scattering  studies  has  not  been  actively  pursued. 

However,  the  relatively  recent  improvement  and  development  of  image  intensifier 
detectors  provides  exactly  the  needed  detection  method.  With  this  type  of  detector,  both 
broad-band  and  high  spectral  resolution  detection  are  possible  for  10-nsec  pulse  sources. 
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This  development  now  provides  the  capability  of  cross  section  measurements,  using  single 
laser  firings,  thus  removing  a  major  obstacle  in  the  use  of  the  shock  tube  for  cross  section 
measurements. 

Having  discussed  the  general  rationale  for  the  use  of  the  shock  tube  for  cross  section 
measurements  and  having  shown  that  recent  instrumentation  developments  have  eliminated 
the  major  objections  to  the  use  of  this  tool,  it  is  necessary  to  consider  both  the  design  criteria 
for  such  a  shock  tube  system  and  to  investigate  the  known  problems  of  operation  of  a 
practical  shock  tube  system.  Subsequent  sections  of  this  report  will  specify  the  operational 
requirements  and  design  criteria  of  the  shock  tube  and  the  considerations  of  first-order 
corrections  for  real-gas  effects,  such  as  shock  wave  attenuation,  boundary-layer  formation, 
and  degradation  of  the  ideal  run  time.  Additionally,  a  description  will  be  presented  of  the 
shock  tube  system  which  was  fabricated  to  satisfy  these  requirements.  Finally,  the  laser 
scattering  system  which  was  assembled  to  perform  the  Raman  scattering  measurements  will 
be  described,  and  preliminary  results  of  vibrational  Raman  scattering  measurements  of 
shock  heated  N2  will  be  presented. 

2.0  DESIGN  CONSIDERATIONS 
2.1  IDEAL  SHOCK  TUBE  PERFORMANCE  PREDICTIONS 

As  a  brief  review  of  the  operation  of  the  shock  tube,  Fig.  2  sketches  the  trajectory  of  an 
ideal  shock  wave  and  defines  the  various  regions  of  the  gas  shock  flow  field.  As  seen  in  the 
depiction  of  the  pressure  wave,  the  parameter  subscripts  denote  the  following  regions  of  the 
shock  tube: 


Subscripts 

Region 

l 

Driven  gas 

2 

Incident  shock  wave  region  be¬ 
tween  shock  front  and  contact 
surface 

3 

Region  between  contact  sur¬ 
face  and  rarefaction  wave 

4 

Dr  iv«r  gas 

5 

Reflected  shock  region 
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The  quantitative  description  of  the  shock  wave  is  given  by  the  following  relations  (Ref. 

12). 


Incident  shock  pressure  ratio: 


P2  =  [2TM?  -  (7-I)] 
Pi  (7  +  l) 


Incident  shock  temperature  ratio 
t2 


12  _  r  m2 

t7  -  PM|-—  Jx 


[ 


— ; — Ml  +  1 


[^M,]: 


Incident  shock  mass  and  number  density  ratios: 

62  _  »2  _  (7+ 1)  M) 

ei  ni  [(7-1)  M? -2] 


Incident  shock  Mach  number: 


Mi 


Ws 

ai 


Reflected  shock  pressure  ratio: 

P5  _  f  [27M?-(y- 1)]  ^  f  [O7-  l)Mi  —  2(7—  1)]  1 

P|  \  (7+1)  /  \  [(7-1)Mi+2]  / 

Reflected  shock  temperature  ratio: 

=  {2(7-  1)M?  +  (3-7)}  {(37-QM?  -  2(7-1)} 

T|  [(7  +  DM1]2 


(1) 


(2) 

(3) 

(4) 

(5) 

(6) 


P5  _  ns 
Pi  ni 


(7) 


Figure  3  shows  the  variation  of  P2/Pi,  T2/Ti,  n2/ni  with  incident  Mach  number  Mi  and 
initial  pressure  ratio  P4/P1.  Figure  4  shows  a  similar  variation  for  T5/T1  and  nj/ni  with  the 
ratio  P4/P1.  For  these  calculations  it  has  been  assumed  that  N2  and  He  are  the  driven  and 
driver  gases,  respectively,  and  that  the  inital  temperatures  of  the  driven  and  driver  gases  are 
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identical.  Further,  for  N2  it  has  been  assumed  that  the  specific  heat  ratio  7  =  1.40;  i.e.,  real- 
gas  effects  are  neglected.  From  Fig.  3  it  is  seen  that  for  an  initial  pressure  ratio  P4/P1  ±  200 
one  achieves  an  incident  shock  Mach  number  Mi  ±  4.1,  temperature  ratio  T2/T1  ±  4.2  and 
a  density  ratio  n2/ni  ±  4.6.  Therefore,  for  T|  =  300  K  and  nj  =  0.286  amagat,  the  incident 
shock  wave  values  are  T2  -  1260  K  and  n2  =  1.3  amagat.  Obviously,  smaller  ratios  of  P4/P1 
produce  correspondingly  lower  values  of  T2  and  n2.  If  the  reflected  shock  wave  is  observed, 
from  Fig.  4  it  is  seen  that,  for  the  above  case,  T5  =  2300  K  and  n5  =  4.1  amagat. 
Consequently,  from  the  results  of  these  simple  computations  and  the  results  of  Figs.  3  and  4, 
one  sees  that  it  is  possible  to  produce  a  wide  range  of  values  of  n  and  T  by  varying  the  ratio 
P4/P1  and  consequently,  the  incident  shock  Mach  number  Mj. 


For  Raman  scattering  cross  section  measurements,  it  was  desired  to  produce  an  LTE  gas 
sample  for  which,  in  the  shocked  gas  region,  n<  lO^cc-1  and  500  <  T  <  3000  K. 
Further,  the  pressure  P4  required  to  produce  the  shock  waves  was  desired  to  be  sufficiently 
low  in  magnitude  that  conventional  laboratory  shock  tube  design  and  techniques  could  be 
employed.  Figures  3  and  4  and  Ref.  12  demonstrate  that  these  design  criteria  are  easily 
satisfied.  An  additional  consideration  for  the  design  criteria  for  the  shock  tube  system  was 
the  run  time  or  data  acquisition  time  available  for  the  measurement.  For  incident  shock 
wave  studies,  run  times  of  10-50  /xsec  were  considered  to  be  adequate  provided  that  a  laser 
source  of  pulse  width  tl<1  h sec  was  used.  Obviously,  a  variety  of  lasers  satisfies  this 
requirement,  including  the  Q-switched  ruby  laser,  the  frequency-doubled  Nd:YAG  laser, 
and  the  excimer  laser  sources.  Further,  for  reflected  shock  wave  studies,  the  available  data 
acquisition  time  is  increased,  and  it  was  desired  that  reflected  shock  run  times  on  the  order 
of  500  fisec  or  greater  exist.  If  this  requirement  is  satisfied,  not  only  can  Q-switched  and 
other  nominal  20  /isec  sources  be  used,  but  also,  if  the  background  radiation  and  laser- 
induced  fluorescence  or  scattering  are  sufficiently  small,  it  is  possible  to  employ  the  more 
energetic  conventional  pulsed  ruby  laser  of  nominal  600  /*sec  laser  pulse  width;  the  use  of 
such  a  source  will  increase  significantly  the  S/N  ratio  of  the  detected  signal. 

Although  of  lesser  importance  for  cross  section  studies  than  the  reflected  shock  mode  of 
operation,  consideration  will  be  given  initially  to  the  less  complex  incident  shock  wave.  For 
this  mode  of  operation,  the  duration  of  the  shocked  gas  flow  time  is  defined  to  be  the 
interval  of  time  between  the  passages  of  the  incident  shock  and  contact  surfaces.  Figure  5 
shows  the  x  -t  trajectory  diagram  for  the  incident  shock  (Ref.  12).  From  Fig.  5  it  is  noted 
that  the  incident,  and  assumed  constant,  shock  speed  is  Ws  and  that  the  contact  surface  is 
moving  with  a  constant  speed  V2.  For  an  arbitrary  location  X  in  the  driven  section,  the  run 
time  At,  of  the  incident  shock  wave  is  shown  in  Fig.  5.  It  is  easily  shown  that  A t\  is  given  by 
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Ati  — 


X 

V2 


X 

Ws 


X.  f  K*yi  —  0  Ml  +  2]  1 
a'  i  [2M,(Mi  - 1)]  ) 


(8a) 

(8b) 


Figure  6  shows  the  variation  of  the  product  At;  ai/X  with  incident  Mach  number  M]  for 
—  1 .40.  If  N2  is  assumed  to  be  the  driven  gas  and  X  is  selected  to  be  3  m,  one  finds  At*  to  be 
approximately  990  /tsec. 


Increasing  the  position  of  observation  increases  At;,  as  Eq.  (8b)  shows,  until  one  reaches  the 
position  XCf  shown  in  Fig.  5,  where  the  reflected  rarefaction  head  intersects  the  more  slowly 
moving  contact  surface.  Clearly,  At,  is  maximum  at  x  =  Xc.  The  coordinates  Xc  and  tc  at 
which  this  intersection  occurs  are  given  by 


r  2ai(M|  - 1)  1 
l  [(yj  +  DMiJ  / 


(9a) 


and 


V  /  V  (7)  +  1)  V  a4  * 


(Mf  -  1) 

Mi 


“  (t4  +  0 
1 2(74  -  1) 


(9b) 


where  X4  is  the  length  of  the  driver  section  which  is  charged  with  a  gas  species  characterized 
by  the  isentropic  sound  speed  04,  and  it  has  been  assumed  that  the  flow  behind  the  contact 
surface  is  supersonic. 


At  x  =  Xc,  the  maximum  run  time  At+  is  given  by 

Art  =  rc  -  ts 


=  Xc  [  —  -  —  ], 

1  v2  ws  J 


(10a) 

(10b) 


where 


Ws  -  Miai 


and 


2  arfMi  -  1) 

f(7,  +  1)  Mj] 
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It  can  be  shown  that 


Art  =  tl. 


2(M2i  -  1)  1 
(7,  +  1)M*  J 


01) 


Assuming  that  the  driver  and  driven  gas  temperatures  are  equal  and  that  74  =  5/3  and  7!  = 
7/5,  evaluation  of  Eqs.  (9)  -  (1 1)  for  X4  =  2.4  m  and  Mi  =  4  yields  Ar^  1.25  msec,  which 
is  consistent  with  the  desired  performance  of  the  tube. 


Considering  the  reflected  shock  region',  the  observation  time  (Arr)  can  be  shown  (Ref.  12) 
to  be  given  by 


Arr  = 


X2 

V2 


Xi 

Ws 


(12) 


where  Xi  is  the  driven  length  and  X2  is  the  distance  from  the  diaphragm  to  the  position  of 
interaction  of  the  reflected  shock  and  the  contact  surface.  Further, 


x  ..  [(A)  -  (*)] 
'[(*)-(*)] 

where  WR  is  the  speed  of  the  reflected  shock  and  is  given  by 


(13) 


-  -  { K„ .  _  (_^)} 

Figure  7  shows  the  variation  of  Arr  with  the  length  X4  and  Xj  of  the  driver  and  driven 
sections,  respectively,  for  an  incident  Mach  3  shock  wave.  Selecting  X4  and  X!  to  be  2.1  m 
and  4.6  m,  respectively,  it  is  seen  that  Arr  =  635  /*sec  is  obtained  as  the  observation  time;  this 
value  of  Arr  corresponds  to  the  nominal  conventional  model  pulse  width  of  the  Holobeam 
Model  620  ruby  laser  which  is  to  be  used  for  the  cross-section  studies.  Using  these  selected 
values  of  X4  and  Xi,  the  variation  of  Arr  with  Mi  was  calculated,  and  the  results  are  shown 
in  Fig.  8.  It  is  seen  from  Fig.  8  that,  as  known  from  the  outset,  the  use  of  a  nominal  20-nsec 
laser  pulse  width  is  adequate  for  reflected  shock  studies,  but,  additionally,  it  is  seen  that, 
subject  to  restrictions  which  were  mentioned  previously  about  background  radiation  and 
other  noise  sources,  the  more  energetic  conventional  mode  ruby  laser  is  capable  of 
utilization  for  reflected  shock  studies  for  the  Mach  number  range  2  s  Mi  <;  3.  Consequently, 
from  these  results  it  is  seen  that  the  desired  observation  times,  as  well  as  the  desired 
temperatures  and  densities,  can  be  achieved  with  a  shock  tube  of  reasonable  practical 
lengths  and  P4/P1  pressure  ratios. 
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2.2  NON-IDEAL  EFFECTS  OF  SHOCK  TUBE  OPERATION 

A  variety  of  effects  degrades  the  performance  of  a  practical  shock  tube  relative  to  the 
behavior  of  the  ideal  system.  Although  the  quantitative  evaluation  of  the  majority  of  these 
effects  must  be  obtained  by  calibration  measurements  of  the  shock  tube,  it  is  informative  to 
list  and  discuss  briefly  some  of  these  effects  and  note  their  relative  importance. 

2.2.1  Boundary-Layer  Formation 

The  interaction  of  the  shock  tube  wall  with  a  real  gas  results  in  the  formation  of  a 
boundary  layer  whose  effect  and  importance  increases  with  decreasing  unit  Reynolds 
number  of  the  gas.  Figure  9  sketches  the  boundary  layer  which  is  formed  initially  at  the 
shock  front  and  continues  its  growth  through  the  region  of  the  shocked  gas.  The  importance 
of  this  phenomenon  can  be  attributed  to  three  major  effects  its  existence  has  on  shock  tube 
studies:  (1)  reduction  in  shock  flow  and  observation  time;  (2)  modification  of  the  conversion 
of  laboratory  to  particle-time  scale;  and  (3)  nonuniform  shock  wave  properties.  Each  of 
these  effects  will  be  considered. 

First,  regarding  the  reduction  in  shock  flow  time,  it  has  been  a  common  observation 
(Ref.  12)  for  shock  tube  studies  that  the  observation  time  of  shock  tube  experiments  was 
significantly  less  than  the  ideal  shock  tube  predictions.  Mirels  (Ref.  14)  has  summarized 
these  findings  and  presented  a  quantitative  explanation  of  the  phenomenon.  The  trajectory 
diagram  of  Fig.  10  shows  that  the  distance  between  the  shock  front  and  the  contact  surface 
increases  for  an  ideal  shock  wave.  However,  as  Mirels  shows,  the  formation  of  the  boundary 
layer  results  in  a  maximum  shock  length  (fm)  as  the  contact  surface  accelerates  to  conserve 
the  mass  of  the  “2”  region  gas  which  is  processed  by  either  the  boundary  layer  or  the  contact 
surface  regions.  Upon  achieving  this  steady-state  condition,  the  shock  wave  proceeds  with  a 
constant  length  (0  of  shocked  gas  for  which  f  =  fmax.  This  non-ideal  behavior  is  sketched  in 
Fig.  10,  and,  using  the  previously-defined  observation  time  parameter  At,  Fig.  10  shows 
clearly,  albeit  qualitatively,  the  reduction  in  At  for  a  real  shock  tube.  Equally  clear,  but 
qualitative,  is  the  dependence  of  the  resulting  test  time  reduction  on  the  degree  of  turbulence 
of  the  boundary  condition.  Mirels  presented  results  based  upon  the  assumptions  of 
boundary  layers  which  were  either  laminar  or  turbulent,  but  no  predictions  were  presented 
for  the  case  for  a  boundary  layer  which  was  characterized  by  the  transition  from  laminar  to 
turbulent.  More  recently.  Bander  and  Sanzone  (Ref.  15)  have  considered  the  problem  and 
have  shown  that  for  such  a  transitional  case  At  can  be  expressed  as 

At  =  ^  +  (Vi)  [(Art)2  +  4  tfi(Arf  -  Art)}.  (15) 
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where  At,  and  A t,-  denote  the  turbulent  and  laminar  test  times,  respectively,  and  trt  is  the 
laboratory  time  at  which  transition  to  turbulence  occurs.  Note  that  the  relation  t^  =  0 
defines  a  fully  turbulent  boundary  layer  and,  for  this  case,  At-  Atx.  For  a  completely 
laminar  flow,  t&  =  At(  =  Ar.  Obviously,  the  crucial  parameter  of  Eq.  (15)  is  t*.  For  the 
shock  wave  conditions  studied  in  Ref.  (15),  it  was  found  that  tft  was  given  by  equating  the 
transition  Reynolds  number  (Re),  and  flow  Reynolds  number  (Re)s 

(Re)t  =  (Re),  (16) 


where  (Re),  was  defined  to  be 


(Re),  = 


V2  tfl  Q2 
*2 


(17) 


where  (Re)s  was 

(Re)s  =  (18) 

^2 

The  shock  tube  diameter  is  D,  and  the  gas  speed,  mass  density,  and  viscosity  are  denoted  by 
V,  q,  and  r\ ,  respectively.  From  Eq.  (16)  and  the  relations  of  Eqs.  (17)  and  (18),  one  finds  the 
nondimensional  time  to  transition,  tft,  to  be  given  by 

7  tft  <l/2)(7,  +  l)Mj 

*  (M?  -  1)  09) 

which  is  the  time  required  for  the  shocked  gas  to  travel  a  length  of  one  tube  diameter. 
Therefore,  if  the  results  of  the  rather  restricted  study  of  Ref.  15  are  valid,  the  test  time  of  a 
shock  tube  flow  is  predictable. 


The  second  major  effect  of  the  boundary  layer  is  the  required  modification  to  the 
magnitude  of  the  particle  time  tp.  The  laboratory  time  tf  of  a  shock  wave  is  the  time  duration 
at  an  observation  location,  as  measured  by  conventional  instrumentation,  between  the 
passage  of  the  shock  front  and  the  time  of  measurement  of  a  property  of  the  gas  sample 
within  the  observation  or  measurement  volume,  i.e.. 


f  d£_ 

0  WS 


f 

WS 


(20) 


In  shock-fixed  coordinates  the  time  (tp)  required  for  a  particle  to  pass  from  the  shock  front 
to  the  distance  (  downstream  of  the  shock  is 


tp 


1  *L 

°  U2 


(21) 


14 


AEDC-TR-83  33 


where  uj  is  the  speed  of  the  gas  in  the  shock-fixed  coordinate  system.  If  uniform  shock 
conditions  exist,  it  is  easily  seen  that 


(22) 


which  is  signficantly  greater  than  tf. 

Since  it  is  the  particle  time  (tp)  which  is  relevant  for  rate  kinetics  and  nonequilibrium 
calculations  for  shock  waves,  it  is  quite  clear  that  knowledge  of  tp  is  essential  for  the  design 
of  experiments  which  are  based  upon  the  assumption  of  LTE  and  yet  which  utilize  gaseous 
molecular  species  for  which  energy  transfer  relaxation  processes  exist.  For  such  studies, 
therefore,  knowledge  of  tp,  supplemented  by  appropriate  rate  kinetics  calculations,  will 
enable  prediction  of  the  required  delay  in  measurement  following  passage  of  the  shock  Front 
to  ensure  LTE  for  the  molecular  species. 

Hobson  and  his  co-workers  (Refs.  16,  17,  and  18)  recognized  the  implications  of  Mirels’ 
results  for  the  calculation  of  t„  and  showed  that  for  a  laminar  boundary  layer 


where 


(23) 


tp/tp  (0)  is  shown  in  Fig.  11  as  a  function  of  £/fm,  and  it  is  seen  that  the  correction  is 
significant.  The  relations  used  for  the  calculation  of  C  are  given  in  Ref.  14  and  will  not  be 
reproduced  here. 

In  addition  to  the  two  boundary-layer  effects  which  have  been  discussed  previously, 
there  exists  a  third  effect  which  can  introduce  systematic  errors  in  cross  section 
measurements  obtained  with  the  shock  tube.  This  third  effect  is  the  existence  of  nonuniform 
gas  properties  between  the  shock  front  and  contact  surface  locations.  It  is  necessary,  as  a 
result  of  this  nonuniformity,  to  be  able  to  predict  values  of  the  shocked  gas  properties  as  a 
function  of  the  gas  sample’s  axial  distance  f  from  the  shock  front.  Mirels  (Ref.  19)  has 
considered  this  complication  and  has  shown  that  for  a  shock  Mach  number  Ms  for  a  gas  of 
specific  heat  ratio  y 

$2 1  _  [  2  +  fry  -  1)  mL  ]  ^ 

*2s  [2  +  (7  +  1)  Ms2  ] 
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where  the  subscripts  “2f”  and“2s”  represent  values  of  <t>  at  f  and  at  the  shock  front, 
respectively.  Further,  $  denotes  T,  nt-1,  and  p<7- ■>/■>.  Further,  as  shown  in  Ref.  19,  the 
Mach  number  of  the  flow  directly  behind  the  shock  in  shock-fixed  coordinates  is 


M2s  = 


/  [(7  -  1)  M?  +  2]  1  1/1 

l  [27M^  -  (7  -  1)]  ) 


(25) 


and  the  Mach  number  M2  at  any  location  f/fm  is  obtained  by  numerical  solution  of  the 
equation 


M2 


['  -  01 x  { 


[2  -  (7  -  1)M:]  'I 

[2  (y  -  1)  Ml  ]  J 


(7  +  1) 

[2(7  -  1)] 


(26) 


To  obtain  an  estimate  of  the  maximum  variation  of  the  gas  parameters  within  the  shock,  use 
is  made  of  the  ratio  given  in  Ref.  19  of  the  properties  at  the  contact  surface  (<b2c)  and 
downstream  of  the  shock  front  (<t2s).  These  results  were  tabulated  by  Mirels  for  y  =  5/3, 
7/5  and  6/5,  and  these  results  are  shown  in  Figs.  12  and  13.  Figure  12  shows  the  temperature 
ratio  T2c/T2s)  and  it  is  seen  that  T2c/T2s  decreases  as  Ms  increases.  For  7  =  5/3,  it  is  seen 
that  the  corrections  for  Ms  =  2  and  Ms  =  50  are  12  percent  and  approximately  7  percent, 
respectively.  For  7  =  7/5,  these  corrections  are  approximately  7  percent  and  4  percent, 
respectively.  As  Eq.  (24)  and  Fig.  13  show,  even  larger  variations  exist  for  the  ratio  n2c/n2s. 
As  an  example,  from  Fig.  13,  n2c/n2s  for  7  =  7/5  are  approximately  18  percent  and  8  percent 
for  Ms  =  2  and  50,  respectively.  Although  <F2/$2s  will  be  less  than  4>2C/4)2S,  it  is  clear  that 
accurate  cross  section  measurements  require  consideration  of  the  effects  of  nonuniform 
shocked  gas  properties,  and,  although  not  presented  in  this  report,  these  results  have  been 
extended  to  predict  the  temporal  variation  of  the  shocked  gas  properties  in  the  reflected 
shock  region. 


2.2.2  Real-Gas  Effects 


The  equations  and  relations  of  the  previous  sections  have  assumed  the  shocked  gas 
species  to  be  ideal,  since  the  specific  heat  ratio  7  and  the  chemical  constituency  were 
assumed  to  be  constant.  However,  as  the  gas  temperature  increases,  the  energy  of  excitation 
of  the  internal  and  electronic  modes  of  motion  of  the  molecule  increases.  Consequently,  7 
decreases  as  T  increases,  which  introduces  real-gas  effects  into  the  derivation  of  the  shock 
tube  relations.  The  qualitative  result  of  these  real-gas  effects  produces  post-shock  values  of 
the  temperature  and  density  which  are  lower  and  higher,  respectively,  than  the  ideal  gas 
results.  Similar  effects  result  from  consideration  of  dissociative  reactions  of  the  shocked  gas 
species. 
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The  excitation  and  dissociative  reaction  processes  described  in  the  previous  paragraph 
occur  with  finite  characteristic  times  or  rate  kinetics  coefficients,  and  recognition  that  these 
processes  do  not  occur  with  infinite  speed  introduces  an  additional  complication  in 
prediction  of  the  shocked  gas  parameters.  The  qualitative  features  of  these  processes  are 
most  easily  appreciated  by  considering  the  vibrational  relaxation  mechanism  for  excitation 
of  the  vibrational  mode  of  a  simple  diatomic  species  such  as  N2.  If  the  vibrational  mode  is 
assumed  to  be  in  local  thermal  equilibrium  and  to  be  characterized  by  a  vibrational 
temperature  Tv  corresponding  to  a  molar  vibrational  energy  Ev,  it  is  known  that  the  time 
variation  of  ETv  is  given  by  (Ref.  20). 

^  -  (—)  [Ev  -  Eve],  (27) 

dt  7\ 

where  rv  is  the  vibrational  relaxation  time  and  £ve  is  the  molar  vibrational  energy  if  the  mode 
were  in  equilibrium  with  the  gas  translation  temperature.  The  vibrational  temperature  Tv  is 
related  to  £v  by  the  relation 

r  _  R0v 

E\'  —  _  n  , 

[exp( — -)  -ll 

L  “  Tv  )  J  (28) 

where  it  is  assumed  that  the  diatomic  molecule  is  characterized  by  a  simple  harmonic 
oscillator  of  characteristic  vibrational  temperature  0V.  R  is  the  universal  gas  constant. 

As  N2  is  processed  by  the  shock  front,  which  is  several  mean  free  paths  in  thickness,  the 
relaxation  time  rv  is  sufficiently  long  and  the  rotational  mode  relaxation  time  rR  is 
sufficiently  short  that  N2  behaves  as  a  gas  with  an  approximate  7  of  7/5.  Consequently,  the 
temperature  and  density  increase  through  the  shock  front  follow,  in  the  main,  the  path  for  a 
7  =  7/5  gas,  and  this  is  shown  for  the  temperature  increase  in  Fig.  14.  Following  this 
increase  in  both  rotational  and  translational  temperatures  through  the  shock,  vibrational 
excitation  becomes  significant,  and  Ev  attempts  to  relax  toward  the  variable  value  of  the 
translational  temperature  Ttr  which  is  assumed  to  equal  the  rotational  temperature  TR;  i.e., 
Ttr  =  TR  =  T. 

Since  an  increase  in  Ev  is  supplied  from  the  translational  and  rotational  modes,  an 
increase  in  Ev  and,  therefore,  Tv,  results  in  decreasing  values  of  Ttr  and  TR,  as  shown  in  Fig. 
14.  Quite  clearly,  if  measurements  of  the  Raman  spectra  of  the  vibrational  modes  of 
molecular  motion  are  to  be  observed,  the  speed  with  which  vibrational  equilibrium  is 
attained  is  a  crucial  parameter  for  cross  section  measurements. 
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Since  tv  for  most  diatomic  species  is  described  by  the  Landau-Teller  theory  which 
predicts  that  tv  varies  as  exp(T~1/3),  rv  decreases  rapidly  as  T  increases  (Ref.  20). 
Consequently,  more  rapid  vibrational  relaxation  is  expected  for  high  Mach  number  shots 
and  for  the  reflected  shock  region.  However,  prediction  of  this  effect  is  essential  for  accurate 
shock  tube  studies.  Although  they  will  not  be  reported  in  this  work,  such  calculations  have 
been  performed  for  N2  and  binary  mixtures  of  N2  and  02,  and  N2  and  H20,  and  the  effect  of 
the  boundary-layer  formation  on  this  collisional  phenomenon  has  been  included  in  the 
computation.  For  this  report,  the  measurements  of  the  vibrational  Raman  spectra  have  been 
restricted  to  temporal  regions  in  the  shock  flow  which  are  either  late  in  the  incident  shock  or 
early  in  the  reflected  shock  regions.  Therefore,  even  though  relaxation  calculations  are 
required  for  accurate  predictions,  it  is  ensured  that  significant  vibrational  excitation  will 
exist. 


3.0  SHOCK  TUBE  SYSTEM 
3.1  SHOCK  TUBE  CONFIGURATION 

Based  on  the  experimental  considerations  which  were  either  discussed  or  implied  in  the 
previous  sections,  the  design  criteria  were  determined,  and  the  operating  condition  ranges 
and  shock  tube  dimensions  were  specified.  The  configuration  chosen  was  a  circular  cross- 
section  tube  of  6.35-cm  diam  with  a  driver  section  length  of  2.4  m.  The  driven  section  length 
was  4.57  m.  The  mismatch  in  diameter  of  the  test/driven  section  was  0.025  mm.  However, 
deviations  from  circularity  of  the  tubes  and  flanges  produced  a  circumferentially-random 
maximum  mismatch  of  0.050  mm.  To  satisfy  the  safety  requirements  for  operation,  the 
driven  section  of  the  tube  was  fabricated  from  1018  seamless  steel  tubing  of  2.75-in.  (6.99 
cm)  OD  and  0.120-in.  (0.30  cm)  wall  thickness.  The  driver  section  was  fabricated  from 
2.5-in.  (6.35)  ID  black  steel  pipe  with  a  0. 250-in.  (0.64  cm)  wall  thickness,  which  has  a  design 
and  statically  tested  maximum  charge  pressure  of  1 ,000  psia  (68  atm).  The  shock  tube  was 
mounted  on  a  12-in.  (30.5  cm)  channel  iron  bar  which,  in  turn,  was  attached  to  five 
4-in. -(10. 2  cm)  diam  posts.  The  posts  were  braced  and  bolted  to  the  concrete  floor.  The 
design  of  the  mounting  system  was  such  that  deflection  of  the  shock  tube  base  during  a  shot 
would  be  less  than  or  equal  to  0,130  mm,  and  alignment  was  accomplished  using  a  He-Ne 
laser  and  level.  Figure  15  is  a  sketch  of  the  shock  tube  and  the  mounting  installation. 

Before  initial  operation  of  the  shock  tube,  it  was  evident  that  the  interior  wall  of  the  tube 
had  rusted.  Operation  of  the  shock  tube  in  conjunction  with  the  laser  scattering  system, 
which  is  described  in  a  following  section,  showed  that  the  particulate  content  of  the  shocked 
gas  was  excessively  large  for  Raman  scattering  studies.  To  eliminate  the  rust  problem  and  to 
minimize  particulate  scrubbing  from  the  wall,  annular  liners  of  304  stainless  steel  were 
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inserted  into  both  the  driver  and  driven  sections  of  the  shock  tube.  The  liner  was  of  0. 188- in. 
(0.476  cm)  thickness  and  was  press-fit  into  the  tube,  and  soldering  was  used  at  the  tube 
penetrations  to  achieve  both  vacuum  integrity  and  stationary  positioning  of  the  liner.  The 
final  inner  diameter  of  the  tube  was  2.125  in.  (5.60  cm). 

Diaphragm  rupture  is  achieved  by  the  conventional  pressure  burst  method  and, 
alternatively,  to  achieve  greater  shock  reproducibility,  if  necessary,  by  the  use  of  a  plunger 
device.  The  plunger  is  located  upstream  of  the  diaphragm  and  is  oriented  at  an  angle  of  37 
deg  with  respect  to  the  tube  centerline.  Actuation  of  the  sharp-tipped  plunger  is 
accomplished  by  air  pressure,  and  12  msec  are  required  to  achieve  the  3.8-cm  extension  of 
the  plunger  from  its  retracted  position  to  the  point  where  it  pierces  the  diaphragm.  The 
diaphragm  materials  used  are  Mylar®  or  annealed  aluminum,  depending  on  the  shock 
strength  desired.  For  modest  shock  strengths,  Mj  =  4,  one  or  more  pieces  of  1-5  mil  (0.03 
-1 3mm)-thick  Mylar®  are  used.  Larger  Mach  numbers  and  pressure  ratios  require  the  use  of 
scribed  metal  diaphragms.  With  this  plunger  system,  the  driver  section  can  be  filled  to  a 
pressure  inferior  to  the  burst  pressure,  the  driver  gas  can  be  stabilized,  and  the  diaphragm 
plunger  can  be  activated. 

The  gas  charging  system  is  divided  into  two  parts  for  filling  the  test  gas  and  driver  gas 
sections.  Provision  has  been  made  to  cold  trap  the  condensible  impurities  of  the  test  gas.  The 
high-pressure  side  is  charged  with  He  for  shock  operations.  Because  the  low-pressure  system 
needs  protection  from  high  pressure,  interlocks  have  been  included  to  inhibit  pressure 
release  which  would  be  hazardous  to  either  personnel  or  equipment.  Firing  is  also  inhibited 
until  all  essential  valves  have  been  closed.  All  firing  and  high-pressure  charging  is  done 
remotely  from  a  control  panel  on  which  are  located  the  power  switches  for  vacuum  pumps, 
thermocouple  readouts,  and  vacuum  gauges.  Nitrogen  is  available  to  purge  the  system  and 
maintain  a  positive  pressure  in  the  tube  during  periods  of  infrequent  use. 

The  shock  tube  is  equipped  with  mechanical  pumps  and  an  oil  diffusion  pump  for  fast 
removal  of  post-shock  gases,  and,  using  these  pumps,  a  pressure  of  approximately  10-5  torr 
is  obtainable  in  both  the  driver  and  driven  sections.  The  diffusion  pump  cold  trap  is  cooled 
to  —  30°C  by  a  Vz- ton  refrigerator  to  minimize  oil  backstreaming,  and,  as  a  result,  the  rate 
of  outgassing  is  approximately  2  mTorr/min.  Consequently,  if  the  test  gas  pressure  is  10  torr 
and  10  min  are  required  to  charge  the  gas  systems  and  initiate  the  shock,  the  impurity  level 
can  approach  2  parts  per  thousand. 

Relief  valves  in  both  the  driver  and  driven  sections  are  essential  design  features  for 
protection  of  the  shock  tube  and  operator.  In  the  high-pressure  section,  a  relief  valve, 
adjusted  for  1,050  psia  (71.4  atm),  does  not  allow  the  design  charge  pressure  to  be  exceeded. 
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The  low-pressure  section  relief  valve  is  preset  at  approximately  350  psia  (23.8  atm)  for 
protection  of  the  quartz  windows  in  the  test  cell. 

Stack  venting  from  both  relief  valves,  as  well  as  all  normal  high-pressure  venting,  is  to 
the  outdoors  through  1-in. -(2. 54  cm)-diam  stainless  steel  tubing.  The  venting  system 
capacity  is  sufficient  to  withstand  a  complete  high-pressure  regulator  failure.  Since  toxic 
gases  are  not  to  be  used  initially  and  the  pumped  volume  is  small,  the  vacuum  pumps  have 
been  vented  to  the  ambient  atmosphere.  Modification  can  be  simply  made  if  stack  venting  is 
required. 

3.2  CONVENTIONAL  INSTRUMENTATION 

The  equations  of  the  theoretical  section  show  clearly  that  knowledge  of  the  shocked  gas 
parmeters  requires  measurement  of  Pi,  P4,  T[,  T4,  and  the  incident  shock  speed  W$. 
Further,  for  operation  of  the  tube  in  either  the  simple  reflected  or  the  tailored  mode, 
measurements  of  P5  and  T5  are  desired.  The  measurements  of  P]  and  P4  are  performed  using 
calibrated  Tabor  gauges  which  have  0. 1-percent  accuracy.  For  the  P4  measurement  a  2,000- 
psia  (136  atm)  gauge  is  used,  and,  at  present,  a  5-psia  (0.34  atm)  gauge  is  used  in  the  driven 
section.  Additional  gauges  in  both  sections  are  used  to  provide  a  visual  display  of  the 
pressure.  The  pressure  gauges,  as  well  as  the  other  conventional  instrumentation  used,  are 
listed  in  Table  1. 

A  Kistler  gauge  was  located  in  the  flange  at  the  end  of  the  driven  tube.  Calibrations  were 
performed  to  obtain  full-scale  readout  at  300  psia  (20.4  atm)  dynamic  pressure.  Expected 
accuracy  is  of  the  order  of  5  percent  of  full  scale  with  a  rise  time  of  approximately  10  jisec. 

To  obtain  the  incident  shock  velocity,  a  series  of  platinum  heat  transfer  (HT)  gauges 
were  installed  upstream  of  the  test  section.  Figure  16  illustrates  the  design.  The  HT  gauges 
are  spaced  uniformly  30.5  cm  apart  with  the  last  gauge  1  m  from  the  end  of  the  tube. 
Engelhard  05-X  platinum  paint  was  brushed  across  a  6.35-mm-diam  quartz  substrate  and 
fired  to  produce  a  strip  of  platinum  approximately  0.5-mm  wide.  Resistances  desired  and 
obtained  for  all  the  strips  were  50  ±50.  The  strips  are  oriented  perpendicular  to  the  shock 
flow  for  maximum  time  resolution.  The  shock  velocity  is  obtained  by  measuring  the  time  of 
transit  of  the  shock  front  between  the  HT  gauges,  which  are  separated  by  known  distances. 
For  a  Mach  4  shock,  the  incident  shock  velocity  is  approximately  1.4  mmy>sec.  The  time 
interval  between  successive  HT  gauges,  assuming  constant  velocity,  is  then  about  427/isec. 
The  uncertainty  in  the  time  interval  due  to  the  width  of  the  strips  is  about  0.7/isec,  or 
approximately  0.2  percent. 
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3.3  LASER  SYSTEM  AND  OBSERVATION  CELL 

Since  the  Raman  cross  section  measurements  are  to  be  performed  in  both  the  incident 
and  reflected  modes,  a  Raman  scattering  observation  cell  has  been  attached,  as  shown  in 
Fig.  17,  to  the  end  of  the  driven  section  of  the  shock  tube.  The  cross  section  of  the  interior  of 
the  scattering  cell  is,  of  course,  circular,  and  the  tolerance  of  0.050  mm  was  maintained  for 
non-coincidence  of  the  circular  cross  sections  of  the  scattering  cell  and  the  driven  tube 
section.  The  aluminum  scattering  cell  includes  four  optical  ports  into  each  of  which  are 
positioned  3.49-cm-diam  quartz  windows  of  6.35-mm  thickness. 

As  shown  in  Fig.  17,  the  laser  beam  is  injected  vertically  downward  into  the  observation 
volume,  and  the  incident  beam  is  collected  and  dissipated  by  a  beam  collector/dump 
assembly  which  is  located  opposite  the  injection  optics  port.  The  basic  laser  source  is  a 
Model  620  Holobeam  pulsed  ruby  laser  which  provides  an  output  at  694.3  nm  with  a 
nominal  pulse  width  and  energy  of  650  gsec  and  36  J,  respectively,  in  the  conventional  mode 
of  operation.  Using  a  Pockels  cell  Q-switch,  the  capability  exists  for  operation  with  a  pulse 
width  of  25  nsec  and  a  pulse  energy  of  3  J.  A  single-stage  amplifier  is  available  for  addition 
to  the  oscillator  to  increase  the  laser  beam  energy  to  a  maximum  value  of  approximately  6  J. 
Because  of  the  observation  of  particulates  in  the  shocked  gas  which  results  in  laser-heated 
particulate  incandescence,  as  described  and  demonstrated  in  Ref.  10,  the  use  of  frequency- 
doubling  of  the  laser  beam  shifts  the  Raman  scattering  signal  to  the  blue-region  tail  of  the 
near-blackbody  particulate  incandescence  spectrum.  Consequently,  the  signal-to- 
incadescence  noise  is  increased  by  approximately  two  orders  of  magnitude.  To  accomplish 
this  wavelength  shift,  a  rubidium  dihydrogen  arsenate  (RDA)  frequency-doubling  crystal 
was  used,  and  the  ruby  laser  beam  at  694.3  nm  was  converted  with  approximately  20  percent 
efficiency  to  a  wavelength  of  347.15  nm. 

For  the  initial  measurements,  the  laser  oscillator  was  operated  at  an  approximate  output 
power  of  500  MW  and  the  frequency-doubled  output  power  was  approximately  100  MW, 
The  0.635-cm-diam  beam  at  347.15  nm  was  focused  onto  the  centerline  of  the  shock  tube 
using  a  lens  system  of  an  effective  focal  length  of  7.62  cm.  The  approximately  cylindrical 
focal  volume  dimensions  were  1-mm-diam  and  6-mm  length  which  yields  an  approximate 
power  areal  density  at  the  scattering  volume  of  10,D  W/cm2.  The  collected  scattered 
radiation  was  focused  onto  the  400-/xm-wide  entrance  slit  of  an  0.85-^m  focal  length,  double 
spectrometer  which  was  found  to  have  a  typical  reciprocal  dispersion  of  0.44  nm/mm.  The 
dispersed  radiation  was  passed  through  an  exit  slit  of  400  *im  onto  an  RCA  31034 
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photomultiplier  tube  (PMT)  which  was  operated  at  2,000  V.  Data  acquisition  and 
registration  were  achieved  using  a  gated  integrator  and  CRT  display  and  photographic 
recording.  Available  for  data  acquisition  and  preliminary  data  reduction  of  the  pressure 
gauges,  the  PMT,  laser  beam  energy  monitor,  and  heat  transfer  gauges  was  an  on-line  PDP- 
8E  minicomputer  system.  More  complete  details  of  the  laser  scattering  system  are  to  be 
found  in  Refs.  1-3,  6,  and  7. 


4.0  RESULTS  AND  CONCLUSIONS 


4.1  RESULTS 

A  series  of  measurements  was  performed  to  demonstrate  the  operation  of  the  shock 
tube,  its  associated  instrumentation,  and  the  laser  Raman  scattering  system  for  temperature 
measurements.  For  this  purpose,  He  was  used  as  the  driver  gas  and  N2,  a  diatomic  species 
with  a  well-known  Raman  scattering  profile,  was  selected  as  the  driven  gas.  The  gas  shock 
speed  Ws  was  determined  using  heat-transfer  gauges,  and  the  variation  of  Ws  with  the 
diaphragm  pressure  ratio  P4/P1  was  determined.  Vibrational  Raman  scattering  of  N2  was 
observed  at  two  spectral  locations  of  the  Stokes  vibration-rotation  band,  and  the  vibrational 
temperature  Tv  was  determined  using  a  spectral  bandpass  such  that  the  Tv  measurements 
were  comparatively  insensitive  to  the  value  of  the  rotational  temperature  Tr. 

The  initial  studies  were  performed  with  the  seamless  steel  shock  tube  using  both  Mylar® 
and  scribed  aluminum  diaphragms.  The  Raman  scattering  data  using  the  frequency-doubled 
ruby  laser  exhibited  an  intolerable  level  of  background  radiation  because  of  particulate  Mie 
scattering  and  laser-induced  particulate  incandescence.  Although  a  fraction  of  this 
background  radiation  could  be  attributed  to  particulate  dust  created  at  diaphragm  rupture, 
it  was  obvious  that  the  overwhelming  majority  of  the  signal  resulted  from  the  rust  particles 
which  were  scrubbed  from  the  interior  tube  wall  by  the  shock  front.  The  most  expeditious 
solution  to” this  problem  was  the  insertion  of  a  stainless  steel  tube  liner  into  both  the  driver 
and  driven  sections  of  the  shock  tube;  the  results  reported  in  this  work  were  obtained  using 
this  tube  liner  which  successfully  reduced  the  background  particulate  component  to 
approximately  10  percent  of  the  Raman  signal. 

Demonstration  of  shock  tube  operation  was  achieved  over  the  incident  shock  Mach 
number  range  of  approximately  4  to  5,  and  the  shock  speed  Ws  was  measured  using  arrival 
times  of  the  incident  shock  front  at  the  known  positions  of  the  heat-transfer  gauges  which 
were  described  in  a  previous  section.  For  operational  purposes,  Ws  was  determined  as  a 
function  of  the  diaphragm  pressure  ratio  P4/P],  which  is  a  primary  independent 
experimental  variable.  Figure  18  shows  partial  results  of  the  variation  of  P4/P1  with  Ws,  and 
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it  is  seen  that  deviations  from  the  theoretical  prediction  of  approximately  2  to  3  percent  in 
Ws  result  from  a  selected  value  of  P4/P1.  However,  these  differences  in  Ws  between  the 
theoretical  and  experimental  are  repeatable  to  within  approximately  1  percent  which 
determines  the  repeatability  of  a  specified  condition. 


For  the  shock  speed  measurement,  erratic  performance  was  encountered  for  one  of  the 
heat-transfer  gauges  and  for  the  end-wall  pressure  gauge,  and  data  were  acquired  without 
their  use.  The  result  of  the  absence  of  the  heat-transfer  gauge  was  slightly  worsened  accuracy 
in  the  measurement  of  the  shock  speed  Ws.  Figure  19  shows  sample  data  of  the  time  history 
of  the  locus  of  the  incident  shock  front.  For  a  single  measurement  of  Tv  of  N2  using  Raman 
scattering,  a  series  of  four  shots  were  performed:  two  on-band  Raman  data  shots  at 
wavelengths  Xq  and  Xt  and  two  off-band  background  shots  at  wavelength  XB.  For  a  single 
series  of  such  shock  tube  firings,  from  the  data  of  Fig.  19  and  the  data  of  the  two  additional 
shots  it  was  found  that 


Ws  =  (1.684  ±  0.022)  x  103  m/s 
=  (1.704  ±  0.008)  x  103  m/s 
=  (1.714  ±  0.008)  x  103  m/s 
=  (1.710  ±  0.004)  x  103  m/s 


where  the  square  root  of  the  variance  of  Ws  for  each  shot  is  quoted  as  the  error  in  Wg.  The 
mean  value  of  Ws  for  this  data  set  of  four  shots  is  Ws  =  (1.706  ±  0.014)  x  103  m/sec. 
Consequently,  for  a  single  shot,  the  heat-transfer  gauge  data  yields  a  shock  speed 
imprecision  of  less  than  approximately  1  percent,  and  the  speed  of  a  series  (Ws)  can  be 
obtained  with  an  imprecision  level  of  about  1  percent.  Consequently,  the  temperature  ratio 
Tj/T,,  which  is  computed  from  Ws  and  Mi,  can  be  expected  to  vary  for  a  series  of  such 
shots  by  approximately  1 .5  -  2.0  percent  for  an  ideal  gas  case  where  y  is  constant,  and  T2/T 1 
for  a  single  shot  will  be  uncertain  by  approximately  0.5  to  1.0  percent.  For  this  series  of 
measurements,  the  mean  value,  Mi,  of  the  incident  Mach  number  is  Mj  =  4.96  ±  0.004. 
The  incident  shock  temperature  T2  is  approximately  1600  K. 


The  malfunction  of  the  end-wall  pressure  gauge  had  serious  implications  for  proper 
timing  of  the  laser  pulse.  If  the  laser  scattering  observation  volume  is  located  sufficiently  far 
from  the  shock  tube  end- wall,  timing  of  the  laser  firing  and  determination  of  the  location  of 
the  laser  measurement  within  the  incident  shock  wave  is  straightforward.  Measurement  of 
Ws  and  the  time  of  laser  firing  with  respect  to  passage  of  the  shock  front  past  any  one  of  the 
HT  gauges  will  be  sufficient  for  an  accurate  location  of  the  measurement  in  the  shock. 
However,  if  the  end-wall  is  not  located  at  a  large  distance  from  the  observation  volume  and 
measurements  are  being  made  in  the  later  regions  of  the  incident  shock,  it  is  essential  to 
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know  when  the  incident  shock  reaches  the  end-wall.  Then,  assuming  an  ideal  gas,  the 
reflected  shock  speed  WR,  which  is  obtained  from  (Ref.  12): 

WS  J  7+1  P2  T 

\  (7-1)  Pi  / 

is  used  to  find  the  time  of  arrival  of  the  reflected  shock  at  the  observation  point.  For  N2,  for 
which  7  varies  from  1.4  in  the  incident  shock  to  less  than  1.3  in  the  reflected  shock,  one 
finds,  for  P2/P,  =  30,  that  V  (WR/WS)  equals  0.363  and  0.291  for  7  =  1.4  and  1.3, 
respectively.  The  timing  of  the  laser  firing  was  such  that  the  later  region  of  the  incident 
shock  was  observed  if  the  reflected  shock  were  behaving  as  a  7  =  1.3  molecule.  However, 
the  20-percent  difference  in  WR/WS  is  sufficiently  large  that,  if  the  vibrationally  relaxing  N2 
is  behaving  at  the  time  of  the  laser  firing  as  a  7  =  1 .4  species,  the  early  region  of  the  reflected 
shock  was  viewed  rather  than  the  end  of  the  incident  shock.  This  uncertainty  would  have 
been  removed  with  accurate  end-wall  timing,  measurements  of  P5,  and  observation  at  the 
test  volume  of  time  of  passage  of  the  reflected  shock  front. 

Vibrational  Raman  scattering  measurements  were  performed  using  a  calibrated  single 
PMT  channel  in  conjunction  with  the  apparatus  described  in  a  previous  section.  To 
determine  the  wavelength  or  spectral  locations  of  these  measurements,  Raman  spectral 
profiles  were  calculated  over  a  range  of  temperatures,  and  Fig.  20  shows  the  results  of  this 
calculation  for  the  selected  temperature  of  2250  K  for  a  0.85-m  focal  length  double 
spectrometer  which  was  operated  with  0.4-mm  entrance  and  exit  slit  widths,  a  1200 
groove/ mm  grating  and  with  a  reciprocal  linear  dispersion  of  0,44  nm/mm.  For  the  range  of 
temperatures  anticipated,  measurements  were  performed  at  the  wavelenths  of  X<,  =  377.66 
nm  and  X,  =  377.26  nm,  corresponding  approximately  to  the  peaks  of  Raman  scattering 
from  the  v  =  0  and  v  =  1  N2  vibrational  energy  levels,  respectively.  The  scattering  signal 
levels  at  Xq  and  Xh  S(X0),  and  S  (X,),  respectively,  are  shown  in  Fig.  21  as  a  function  of  the 
temperature,  and,  for  the  results  shown  in  both  Figs.  20  and  21,  it  was  assumed  that  the  N2 
gas  sample  exhibited  local  thermal  equilibrium;  i.e.,  T  =  Tv  was  assumed.  Evident  in  Fig.  21 
is  the  depletion  in  number  density  population  of  the  v  =  0  level  and  the  increase  in 
population  of  the  v  =  1  level  as  the  temperature  increases.  Figure  22  shows  the  ratio  Rv  = 
S(Xj)/S(Xo)  as  a  function  of  temperature,  and  the  approximately  linear  dependence  of  Rv  on 
Tv  is  to  be  noted. 

Raman  scattering  data  were  acquired  at  the  wavelength  calibrated  spectral  locations  of 
and  Xj.  Successive  shock  tube  firings  were  required  for  these  measurements,  and  a  series  of 
shots  was  performed  to  obtain  data  at  X0,  Xj ,  and  the  background  radiation  location  Xg  = 
378.3nm  where  only  background  and  particulate  radiation  components  existed.  These 
background  data  indicated  signals  of  approximately  10  percent  of  the  magnitude  of  the 
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Raman  signals,  and  simple  subtraction  was  used  to  correct  the  Raman  signals  for  this 
background  continuum  which  was  assumed  to  vary  little  over  the  spectral  band  pass  of  the 
Raman  measurements.  Figure  23  shows  typical  oscilloscope  traces  for  the  Raman  signals  at 
and  \|,  and  the  nominal  20  nsec  width  and  signal  levels  exceeding  20  mv  are  to  be  noted. 
For  a  series  of  measurements  corresponding  to  Mt  =  4.96  +  0.04,  the  values  of  Tv  were 
determined  and  the  results  were  Tv  =  21 10  ±  220  K,  where  the  uncertainty  value  of  220  K 
represents  the  standard  deviation  of  the  sample. 

The  measured  value  of  Tv  exceeds  the  gas  temperature  T2  =  1600  K  which  was 
determined  by  Ws  measurements,  and  the  approximate  20  percent  difference  exceeds  the 
estimated  uncertainty  in  both  Tv  and  T2,  the  latter  of  which  is  on  the  order  of  2  to  3  percent. 
The  most  likely  cause  of  the  discrepancy  between  Tv  and  T2  is  that  measurements  of  Tv  were 
performed  early  in  the  reflected  shock  region  where  Tv  is  increasing  to  its  final  value,  which 
is  significantly  larger  than  both  1600  K  and  2110  K.  Finally,  the  approximately  10  percent 
error  in  Tv  is  reasonably  consistent  with  the  overall  uncertainty  of  about  7.5  percent  for  S(A<,) 
and  S\]). 

4.2  CONCLUSIONS 

The  demonstration  of  the  laser  Raman-shock  tube  system  has  been  accomplished,  and 
the  ability  to  perform  Raman  measurements  in  a  shock  tube  environment  has  been  shown. 
The  necessity  for  operation  with  a  stainless  steel  interior  wall,  rather  than  seamless  steel,  was 
shown  conclusively,  as  was  the  comparatively  low  levels  of  particulate-related  scattering 
signals  with  the  stainless  steel  surface.  The  vibrational  temperature  Tv  was  found  to  be  21 10 
±  220  K  for  a  series  of  measurements.  These  results  were  acquired  using  a  single  PMT 
channel  for  detection,  and  successive  laser-shock  tube  shots  were  required  to  obtain  the  two 
signals  required  for  the  Tv  determining  ratio.  One  or  more  background  shots  were  necessary 
to  correct  for  the  background,  radiation  effects.  The  uncertainty  in  Tv  could  be  lessened  by 
performing  at  least  two  of  these  measurements  during  a  single  firing  of  the  laser  and  shock 
tube.  It  is  believed  that  uncertainty  levels  of  ±  5  percent  are  obtainable  in  this  manner. 

The  incident  shock  speed  was  determined  using  heat-transfer  gauges,  and  it  was  found 
that  the  shock  speed  Ws  could  be  measured  with  an  imprecision  of  less  than  approximately  1 
percent  and  that  shot-to-shot  repeatability,  as  reflected  by  Ws  data,  was  on  the  order  of  1  to 
2  percent.  Improvements  in  diaphragm  rupture  techniques  could  lessen  this  uncertainty.  The 
shock  Mach  number  Mi  was  found  to  be  4.96  ±  0.04  for  the  series  of  Tv  measurements  and 
this  value  of  M|  corresponds  to  an  ideal  gas  temperature  of  approximately  1600  K.  The 
necessity  for  accurate  diagnostics  of  the  reflected  shock  pressure  and  speed  was  made  clear 
during  the  course  of  this  study.  Failure  to  obtain  such  data  resulted  in  an  uncertainty  in  the 
location  within  the  shock  of  the  laser  scattering  data.  The  discrepancy  between  T2  =  1600  K 
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and  Tv  =  21 10  K  is  suggested  to  arise  from  the  acquisition  of  laser  data  to  have  occurred  in 
the  higher  temperature  reflected  shock  region  rather  than  the  incident  shock  region. 

Use  of  the  laser-shock  tube  system  for  accurate  cross  section  and  calibration 
measurements  will  require  making  the  relatively  minor  and  previously-mentioned 
experimental  improvements  in  the  system  and  incorporating  into  the  theoretical  analysis 
procedure  the  computational  technique  which  predicts  the  effect  of  the  relaxation  processes 
and  the  boundary  layer  on  the  gas  properties.  These  improvements  are  in  progress  at  this 
time. 
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Diaphragm  Test  Cell 

Figure  1.  Typical  shock  tube  configuration  for  reflected  mode  operation. 
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Figure  2.  Diagram  of  shock  tube  regions  at  time  t  after  diaphragm  rupture. 
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Figure  3.  Ideal  gas  temperature,  pressure,  and  number  density  ratios  for  nitrogren  shocks. 


AEDC-TR-83-33 


AEDC-TR-83-33 


Diaphragm  Distance,  x 

Location 

Figure  5.  Shock  tube  time  *  position  loci  diagram 
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Figure  6.  Incident  shock  flow  time  for  N2  versus 
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Figure  15.  Shock  tube  and  mounting  installation 
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Figure  18.  Pressure  ratio  P4/P1  versus  shock  speed  W$. 
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Figure  20.  Raman  spectral  profile  for  Nj. 
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Figure  23.  Oscilloscope  traces  of  Raman  signals  S(X€)  and  S(Xt). 
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Table  1.  Conventional  Instrumentation 


Type 

Nominal  Range 

Parameter  Measured 

Pressure  Transducers 

Tabor 

Tabor 

Kistler 

<2,000  psia  ( 136  atm) 

<5  psia  (0.34  atm) 

<300  psia  (.20  atm) 

Driver  Gas  Pressure 

Test  Gas  Pressure 

Shocked  Gas  Pressure 

Low-Pressure  Gauges 

Thermal  Conductivity  (2) 
Wallace  and  Tiernan 

Ion  Gauge 

^  10-2  torr 

<  10-3  torr 

<  100 torr 

Pumping  Line  Pressures 

Test  Section  Pressure 

Same 

Thermocouples 

Copper  -  Constantan  i2) 

=  250  - 300  K 

Driver  and  Driven  Section 
Temperatures 

Thermistor 

=  250  -  300  K 

Test  Section  Temperatures 

High-Pressure  Gauges 
Bourdon 

1  Bourdon 

<2,000  psia  1 136  atm' 

<600  psia  1 4 1  atmi 

Driver  Gas  Pressure 

Driven  Gas  Pressure 
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ai 


D 

dcr/d  ft 

<“'w 

Ev 

h 

f 


Mi 

Mi 


ni 


ot 


Pi 

P0 

R 

Rv 

(Re),,, 

S(X) 

S/N 


T,Tr,Tv 


NOMENCLATURE 
Speed  of  sound  in  region  1 
Diameter  of  shock  tube 
Differential  cross  section 

Molar  vibrational  energy;  Eve  is  the  equilibrium  value 
Planck’s  constant 
Length  parameters 

Maximum  length  of  boundary-layer  convected  shocked  gas 
Mach  number  in  region  i 

Mean  value  of  experimental  measurements  for  incident  Mach  number 

Number  density  in  region  i 

Total  number  density 

Pressure  in  region  i 

Power  of  laser 

Universal  gas  constant 

S(XL)/S(Xo) 

Reynolds  number  of  flow  and  transition,  respectively 
Raman  signal  at  wavelength  \ 

Signal-to-noise  ratio 

Static,  rotational,  and  vibrational  temperatures,  respectively 
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TV 

Tj 

t 

tf 

lp 

t(0) 

p 


tA 


Ui 

Vi 

WR 

WS 


x 

Xi 


Ti 


V 

X 

<J> 

f'o 

Q 


Mean  value  of  experimental  measurements  of  Tv 

Temperature  in  region  i 

Time 

Laboratory  time 
Particle  time 

Ideal  particle  time 

Laminar-turbulent  boundary-layer  transition  time 
Particle  speed  in  region  i  in  shock  wave  coordinates 
Particle  speed  in  region  i  in  laboratory  coordinates 
Shock  speed  in  reflected  shock  region 
Incident  shock  speed 
Position 

Length  parameters 

Specific  heat  ratio  for  gas  in  region  i 

Gas  viscosity 

Wavelength;  \0  and  X],  are  Raman  signal  wavelengths  and  Xg  is  the 
wavelength  for  background  radiation  measurement 

Parameter  representing  T,  n-v-i,  and  pb-D/i 

Laser  frequency 

Gas  mass  density 
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Characteristic  vibrational  temperature  of  molecular  vibrational  mode 
Defined  by  Eq.  (8b) 

Rotational  relaxation  time 
Vibrational  relaxation  time 

Observation  time  of  incident  shock  wave 
Maximum  observation  time  of  incident  shock  wave 

Observation  time  for  laminar  boundary  layer 
Observation  time  in  reflected  shock  region 
Observation  time  for  turbulent  boundary  layer 


Subtended  solid  angle 


